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ABSTRACT 

We show that by appealing to a Quark-Nova (the explosive transition of a neutron star to a quark 
star) occurring in an He-HMXB system we can account for the liglitcurve of the first superluminous 
SN, DES13S2cmm, discovered by the Dark-Energy Survey. The neutron star’s explosive conversion 
is triggered as a result of accretion during the He-HMXB second Common Envelope phase. The 
dense, relativistic, Quark-Nova ejecta in turn energizes the extended He-rich Common Envelope in 
an inside-out shock heating process. We find an excellent fit (reduced y 2 of 1.09) to the bolometric 
light-curve of SN DES13S2cnnn including the late time emission, which we attribute to Black Hole 
accretion following the conversion of the Quark Star to a Black hole. 

Subject headings: circumstellar matter stars: evolution stars: winds, outflows supernovae: general 
supernovae: individual (SN DES13S2cmm) 


1. INTRODUCTION 

DES13S2cmm is the first confirmed superlumi- 
nous SN (SLSN) f rom the Dark Energy Survey 
dPapadopoulos et abl (120151) L Its bolometric light-curve 
(LC) shows a slowly declining tail when compared to 
other SLSNe. It also shows an initial hump (the first 
four data p oi nts pr ior t o the peak; see Figure 1). 
iPapadopoulos et abl (120151) explored two possible power 
sources ( 56 Ni decay and magnetar spin-down power) but 
both models poorly fit the LC. In this paper we show that 
a Quark-Nova (QN) occurring in a Helium-High-Mass 
X-ray binary (He-HMXB), which experiences a second 
common envelope (CE) phase, can account for the LC 
features of SN DES13S2cmm including the tail. The pa¬ 
per is organized as follows: in §2 we give a brief overview 
of the QN model and the occurrence of QNe in binary 
systems. In §3 we show the results of applying the QN in 
a He-HMXB system to SN DES13S2cmm. We provide a 
brief discussion and some predictions in §4. 

2. QUARK NOVA (QN) MODEL 

The QN is the explosive transition of a neutron star 
(NS) to a quark star (QS) (Ouyed et al. 2002). The 
conversion energy combined with the core collapse of the 
parent NS results in the ejection of Mqn ~ 10 -3 M 0 
of neutron-rich material (Keranen et al. 2005; Ouyed 
& Leahy 2009; Niebergal et al. 2010). This relativistic 
ejecta has a kinetic energy, Uqn, exceeding 10 52 erg and 
an average Lorentz factor Tqn ~ 10. The QN can occur 
following the explosion of a massive single star or in a 
binary system following accretion from the companion. 
The key constraint is for the NS to reach deconfinement 
densities in its core by becoming massive enough (e.g. 
via fall-back during the SN explosion or accretion from 
a companion) or via spin-down (Staff et al. 2006). We 
define Mns,c. as the NS critical mass above which quark 
decon fi nemen t occurs in the NS triggering the QN (see 
lOuved et al.l (|2013al ) for a recent review). Mns.c. varies 
from ~ 1.6M 0 up to 2M 0 and higher for different equa¬ 
tions of state (e.g. Staff et al. 2006). However we as¬ 
sume Mns,c = 2M 0 to account for the recent heavy NS 
observed by Demorest et al. (2010); ~ 2Mg (or heav¬ 


ier) quark stars can exist when one considers interacting 
quarks (e.g. Alford et al. 2007). 

2.1. QNe in single-star systems: dual-shock QNe 

A dual-shock QN (dsQN) happens when the QN oc¬ 
curs days to weeks after the initial SN. The delay allows 
the QN ejecta to catch up to and collide with the SN 
ejecta after it has expanded to large radii (Leahy&Ouyed 
2008; Ouyed et al. 2009). The SN provides the mate¬ 
rial at large radius and the QN re-energizes it, causing 
a re-brightening of the SN. For time delays of days or 
less, the radius of the SN ejecta is small enough that 
only a modest re-brightening results when the QN ejecta 
hits the SN ejecta. However, the neutron-rich QN ejecta 
lead to unique nuclear spallation signatures (Ouyed et al. 
2011a). For longer time-delays, the radius and density of 
the SN ejecta are such that extreme re-brightening is ob¬ 
served which could expla i n some SLSNe (lOuved et all 
l2012t iKostka et al.l 12014 lOuved fe Leahvl 120131) . For 
time-delays exceeding many weeks, the SN ejecta will 
be too large and diffuse to experience any re-brightening. 
The dsQN model has been successfully applied to a num¬ 
ber of superluminous and double-humped supernovae 
(see http://www.quarknova.ca/LCGallery.html for a pic¬ 
ture gallery of the fits). 

2.2. Quark-Novae in binary systems 

A QN will more often occur in tight binaries. To reach 
AfNS.c. and experience a QN event, the NS needs to ac¬ 
crete from a Roche Lobe (RL) overflowing companion 
(Ouyed et al. 2011b&c; Ouyed & Staff 2013; see also 
Ouyed et al. 2015) or during a CE phase of the binary 
as described in this work. A QN in a binary provides: 

(i) A means (the relativistic QN ejecta) to shock, com¬ 
press, and heat the NS companion or the CE. In the 
case of a QN occurring inside a CE (which is considered 
in this paper), the resulting envelope thermal energy is 
given by -EcE.th. = Csh.-UQN which yields an initial shock 
temperature Tq N , s1i . ~ Csh. x 10 9 K x Uq N ,52 x (M q /M c e) 
for a He-rich envelope of mass Mqe and mean molec¬ 
ular weight hce = 4/3; -Sqm,52 is the QN kinetic en¬ 
ergy in units of 10 52 ergs while </ s h. is the shock heat- 
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ing efficiency. The CE kinetic energy is found from 
Eqn = Ece,k + EcE.th. where the CE thermal energy 
is EcE.th. = E C E,gas + EcE.rad., shared between the 
gas and the radiation. We find that even for extreme 
shock efficiency Qh. ~ 1 and for high compression ratio 
(~ 4Tqn ~ 40; see Ouyed&Staff 2013), the timescale re- 
qu ired to bu rn He to O (e.g. using burning rates from 
iHuang fe Yul (l!998fl l is too large compare to the adia¬ 
batic expansion timescale of the CE. 

(ii) The spin-down power from the QN compact rem¬ 
nant (the QS), which provides an additional energy 
source. The QS is born with a magnetic field of the order 
of 10 15 G due to color ferromagnetism inherent to quark 
matter during the transition (Iwazaki 2005). In the sce¬ 
nario considered here where the QS accretes during a CE 
phase, the high accretion rates would lead to QS spin-up 
instead of spin-down and the QS would gain mass until 
it becomes a BH without entering a spin-down phase. 
Thus in the picture presented here, spin-down power can 
be ignored. 

(iii) The presence of a gravitational point mass (the 
QS) which may slow down and trap some of the ejecta 
and provide accretion power (e.g. Ouyed et al. 2011b&c; 
Ouyed et al. 2014; Ouyed et al. 2015). The conversion 
of the QS to a BH during the CE phase provides an addi¬ 
tional source of energy (the BH-accretion phase) which 
powers the LC in addition to energization by the QN 
shock. 

2.2.1. The Quark-Nova in a He-HMXB system 

In previous papers we considered the detonation of 
a CO White Dwarf by the QN (what we called a QN- 
Ia; Ouyed&Staff 2013; Ouyed et al. 2014; Ouyed et al. 
2015). In a QN-Ia, the system experiences a runaway 
mass transfer leading to the formation of a CO-rich torus 
surrounding, and in the close vicinity of, the exploding 
NS. In this paper we consider a QN in a He-HMXB. 
Such a binary is expected to form when an O/B-HMXB 
binary evolves through a CE phase of the supergiant stel¬ 
lar component. If the binary survives this first CE event, 
the resultant system would contain a He-rich core (Hall& 
Tout 2014) with a NS in a relatively tight orbit. For a 
large mass ratio 1 q, a runaway mass transfer is expected 
once the He-rich core overflows its RL leading to the on¬ 
set of a second CE phase 2 . The NS accretes from two 
mass reservoirs (the first CE and the second CE). For a 
NS born with an initial mass of ~ 1.6M 0 , accretion of 
~ O.4M 0 during the two CE phases is necessary to drive 
it above Mns.c.; accreting ~ O.2M 0 of material per CE 
phase is not unreasonable and should be enough to eject 
the first CE (jBrownl (|1995H : lArmitage fc Livid|200fD L If 
the in-spiralling NS accretes enough mass during the sec¬ 
ond CE phase, it should drive the NS above above Mns.c. 
to undergo a QN explosion in an extended CE; the CE 
expands radially outward as the NS in-spirals towards 

1 E.g. q > 3.5 or orbital period P or b. <0.1 days (Ivanova et 
al. 2003). Alternate conditions leading to a second CE phase give 
2.4 < q < 2.7 and P or b. < 0.25 days (Dewi & Pols 2003). See Tauris 
et al. (2015) for other outcomes and for a recent investigation of 
the evolution of these systems. 

2 In Ouyed et al. (2015), where the NS companion is a degen¬ 

erate WD, the runaway accretion is due to companion expansion 

(since R oc M ~& 3 ) which leads to high-density (degenerate) torus 

formation. 


the core. Following the QN explosion the QS continues 
to in-spiral into the core while gaining mass. Below we 
show that if the QS turns into a BH, BH-accretion can 
power the slowly declining tail of DES13S2cmm. A com¬ 
bination of QN shock heating and BH-accretion provides 
the best fit to DES13S2cmm. 

3. APPLICATION TO DES13S2CMM 

The free parameters in our model are (see Table 1): 
the CE mass (Me e), its initial (extended) radius I?ce,o 
at the time of QN shock breakout; the CE expansion 
velocity uqn induced by the QN shock and QN shock 
heating per particle (3/2)/cBTQN, s h.; ciQN.sh., the shock 
propagation parameter (see Appendix). There are three 
parameters for the BH-accretion model based on the pre¬ 
scription of lDexterfcKaseU ( 20131 (see their appendix A): 
Lq (erg s^ 1 * ), yo and n which together define the injec¬ 
tion power L(t) = Lo(t/to)~ n with to = yotd and td the 
photon diffusion time in the CE; BH-accretion turns on 
at t = to so that L(t) = 0 for t < to- Table 1 shows the 
model’s best-fit parameters for the BH-accretion model 
alone and the for the two-component QN+BH-accretion 
model. The best-fits reduced y 2 values are given in the 
last column of the table. 

3.1. BH-accretion fit 

This model does not have the QN explosion if the NS 
experiences a transition to a BH directly while the system 
is in the second CE phase. The NS could turn into a BH 
prior to or during merging with the CE core which forms 
an accretion disk around the BH to power the LC. The 
best-fit is obtained for td = 13.7 days, yo = 0.1 (i.e. 
to = 5.6 days), L 0 = 6 x 10 44 erg s _1 and n = 0.7. 
The fit is shown in Figure [l] along-side the observations 
from IPanadonoulos et all ( 2015 ). The resulting best-fit 
reduced Xred i s ~ 2.77. The initial hump is not well fit 
with this model. 

3.2. QN+BH-accretion fit 

Adding the QN shock heating yields a significant im¬ 
provement in the LC fit with a reduced y 2 ed of ~ 1.09 
(see Table 1). The BH-accretion phase is delayed from 
the QN event by the time required for the QS to turn into 
a BH, merge with the core and trigger accretion. This 
time delay is to at which point the CE has extended to a 
radius Rce.o+tqn^o; Rce.o is the CE radius at QN shock 
breakout (see Appendix). The LC fit of DES13S2cmm 
is s hown in Figure [2] along-side the observations from 
iPapadopoulos et al.1 l2015t ). 

The QN shock energizes the He CE (of mass Me e = 
2 Mq and radius Rce.o = 135Oi? 0 ) and yields the ini¬ 
tial bright and short-lived hump; the corresponding ini¬ 
tial CE temperature Tce.o is calculated from the QN 
shock heating (Tqn.sIi. = 3.7 x 10 6 K; i.e. Csh. ~ 10 -2 
for Mce = 2M 0 ) by including radiation energy den¬ 
sity (see Appendix). The photon diffusion timescale 
is td = (2/3) \/McF.KT b. /(O.30 cvqei) = 12.6 days with 
/? = 13.8 (lArnettl (11982D ): c is the speed of light and KTh. 
the Thompson cross-section. 

4. DISCUSSION AND CONCLUSION 

In this paper we showed that the LC of SN 
DES13S2cmm is best fit with the QN+BH-accretion 
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TABLE 1 

Best-fit parameters for the DES13S2cmm LC in our model. 



He-rich (i.e. second) CE 

QN 

BH Accretion 

| Fit 


Mce (Mq) 

Rce,o(Ro) 

VQN (km/s) 

?QN,sh. (K) 

«QN,sh. 

t d (days) 

yo 

Lo (erg/s) 

n 


BH-accretion 

- 

- 

- 

- 

- 

13.7 

0.10 

6 x 10 44 

0.7 

2.77 

QN+BH- accretion 

2 

1350 

40,000 

3.7 x 10 B 

3/2 

t 

0.45 n 

2.8 x 10 44 

0.8 

1.09 


t Here, t<j is given by t r \ = (2/3)-^/McE^th/ (0.3/3cuqn) = 12.6 days. t'This corresponds to to = yot <\ — 5.6 days. 
The kinetic energy of the CE/ejecta after the QN is Ece,k = (I/2)— 3.2 X 10 52 ergs. 



Fig. 1.— The BH-accretion model fit (solid line) to the light 
curve of SN DES13S2c mm. The observations (the s olid circles and 
the error bars) are from IPapadopoulos et all (120151 1 ■ 



Fig. 2.— The QN+BH-accretion model fit (solid line) to the 
light curve of SN DES13S2cmm. 

model. In our model, the QN occurs during the sec¬ 
ond CE phase of He-HMXB system followed by a BH- 
accretion phase after the QS turns into a BH in the core 
of the CE. The QN shock re-energizes the extended CE 
explaining the initial hump while BH-accretion power 
nicely fits the slowly declining tail. No nuclear burning is 
triggered by the QN shock which means that there should 
be little or no He-burning products in SN DES13S2cmm 
and similar explosions. 

The lack of He-HMXB systems when compared to the 
obser ved large po p ulat ion of HMXBs with Be-type stars 
(e.g. iRaeuzova fc Ponovl (|2005l) and reference therein) 


has been used as an observational argument in favor of 
mergers during the dynamically unstable mass transfer 
phase of the progenitors (e.g. Linden et al. 2012; see 
also van den Heuvel 1976). However, even in extreme 
cases, the theoretically predicted fraction of He-HMXBs 
surviving the CE p has e is s till higher than the observed 
number ( Linden et al.ll20T2t ). If the NS is born massive, 
it will likely accrete enough matter during the first CE 
phase to reach Mns.c.- The resulting QN could remove 
enough matter to unbind the system and bypass the pro¬ 
duction of He-HMXBs. We thus speculate that the QN 
may be partly responsible for the rarity of He-HMXBs. 

If the NS does not accrete enough mass to go QN in 
the second CE phase, it may reach the center and form a 
Thorne-Zytkow object dThornefcZ’vtkowl (|1977l) l. How¬ 
ever, continued accretion should trigger a QN leading to 
the same outcome. I.e. Thorne-Zytkow objects could be 
short-lived. On the other hand, extreme accretion rates 
could turn the NS directly into a BH. 

The idea of a QN in binaries has proven success¬ 
ful in accounting for some features of SNe-Ia (Ouyed 
et al. 2014; Ouyed et al. 2015) and Gamma Ray 
Bursters (Ouyed et al. 2011c) and has been able to 
account for the LC of DES13S2cmm as shown here. 
Its ability in fitting properties of unusual SNe (see 
http://www.quarknova.ca/LCGallery.html) suggest that 
QNe may be an integral part of binary evolution; the QN 
could lead to novel and interesting evolutionary paths. 
Nevertheless, as we have stated before, our model re¬ 
lies on the feasibility of the QN explosion which re¬ 
quires sophisticated simulations of the burning of a NS 
to a QS which are being pursued. Preliminary simu¬ 
lations with consistent treatment of nuclear and neu¬ 
trino reactions, diffusion, and hydrodynamics show in¬ 
stabilities that could lead to a detonation (Niebergal 
et al. 2010; see also Herzo g fe Ropke (2011) and 
lAlbarracin ManriauefeLugon 3 (i2015l )L We have al¬ 
ready argued that a “core-collapse” QN could also re¬ 
sult from the collapse of the quark matter core (Ouyed 
et al. 2013a) which provides another avenue for the QN 
explosion. 


We thank the referee for comments that helped im¬ 
prove this paper. We also thank Thomas Tauris for 
valuable discussion. This work is funded by the Natural 
Sciences and Engineering Research Council of Canada. 
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APPENDIX 


Due to the outward diffusion of photons, the photosphere is moving inward in mass coordinates, slowly at first but 
faster as the density decreases in time. The ejecta interior to the photosphere we refer to as the core. We will assume 
that the thermal energy in the exposed mass in the photosphere (as the cooling front creeps inward) is promptly 
radiated. The interplay between uniform expansion and radiation diffusion defines the evolution of the photosphere as 

-Rphot.(i) = RcE(t) — D{t ) , (1) 

where Resit) = Rce,o + fQN t and Rce,o the CE envelope radius at QN shock breakout (which corresponds to t = 0 
in our model). Here D(t ) is the diffusion length 


D{tf = Dl + 


c 

?1'CE(t)0Th. 


( 2 ) 


where ncE(t) = Nce/Vce( t) is the number density in the CE. The total number of particles in the CE is Nqe = 
(Mce/mceWh) while Vce (t) = (47r/3)i?cE(t) 3 is the volume extended by the CE and tor the Hydrogen atomic mass. 
We define Do as the initial diffusion length scale by setting nc.E,o&Th.Do — 1 where ncE,o = -Vce/Vce, o and the initial 
volume Vce,o = ( 47 r/ 3 )i?^ E 0 . 

The initial QN shock heating per particle (3/2)fcBT , QN,sh. is a free parameter; fce is the Boltzmann constant. The 
heat is redistributed between gas and radiation to get the post-shock CE temperature Tce,o- The relevant equation 
is |fcBncE,oTcE,o + arad?c E 0 = §fcB?iCE,oT , QN,sh., with uce,o the number density (of electrons and ions) and a ra d the 
radiation constant. 

The subsequent evolution of the CE core temperature after the CE is fully shocked is given by T cole (t) = 
TcE,o(RcE,o/RcE{t)) 2 ~ aQ ^' t ‘ h ' ■ To account for a non-uniform initial temperature, we introduce aQN,sh. which pa¬ 
rameterizes complex shock physics beyond the scope of this work. aQN,sh. = 0 corresponds to the case of an adiabatic 
expansion with spatially uniform initial Tce,o- I.e. the internal energy includes only gas internal energy with 7 = 5/3, 
so as time increases, T cole (t) oc -Rce(£) - 2 - ctQN.sh. also allows to account for the presence of radiation, e.g. for spatially 
uniform Tce.o and 7 = 4/3, one uses ciQN.sh. = 1- ciQN.sh. > 0 also can correspond to a radially decreasing initial CE 
temperature so that T core decreases more slowly than it would be for uniform Tce.o- 

The corresponding luminosity is 

Tqn(£) = cv,tot.(i)AT core (t)ncE(i)47ri? p hot.(i) 2 —> (3) 

where the total specific heat is cv.tot.(i) = cy,gas + cy.rad. = §&b + "'“ncdt)^ ■ Here, AT core ~ T core since the 
photosphere cools promptly (i.e. cooling time is much less than the diffusion timescale). The inward photosplieric 
velocity in mass coordinates is dD(t)/dt. When ncE(t) is low, the T 3 term dominates so that the temperature is much 
lower than in the pure gas model. As ncE(t) increases, the gas energy density becomes more important and T rises. 





